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Guidance notes 


This computer book contains those sections of the chapters in Block D 
which require you to use Mathcad. Each of these chapters contains 
instructions as to when you should first refer to particular material in this 
computer book, so you are advised not to work on the activities here until 
you have reached the appropriate points in the chapters. 


In order to use this computer book, you will need the following Mathcad 


files. 


Chapter D1 

221D1-01 Cartesian and polar forms 
221D1-02 Roots of unity 

221D1-03 Complex exponentials (Optional) 


Chapter D2 


221D2-01 Remainders of powers 
221D2-02 Modular arithmetic 


Chapter D3 

221D3-01 Piecewise linear paths (Optional) 
221D3-02 PLPs with symmetry (Optional) 
221D3-03 A rose window (Optional) 


Instructions for installing these files onto your computer’s hard disk, and 
for opening them, are given in Chapter AO of MST121. 


Activities based on software vary both in nature and in length. Sometimes 
the instructions for an activity appear only in the computer book; in other 
cases, instructions are given in the computer book and on screen. 
Feedback on an activity is sometimes provided on screen and sometimes 
given in the computer book. 


For advice on how each computer session fits into suggested study 
patterns, refer to the Study guides in the chapters. 


Note that there are no specified computer activities associated with 
Chapter D4. 


Chapter D1, Section 6 
Complex numbers and Mathcad 


6.1 Complex arithmetic in Mathcad 


Certain quadratic equations have roots that are complex numbers, and you 
may have already seen Mathcad produce such roots in this context. For 
example, suppose that you use the symbolic keyword ‘solve’ to solve the 
equation x? — 6x + 25 = 0. Mathcad gives the two solutions as 


3+ 41 
a4) 


This is as we would expect, since the equation has the two solutions 3 + 42. 
Note that Mathcad uses the usual symbol 7 for /—1. 


Mathcad will perform various manipulations with complex numbers. ‘To 
use it for this purpose, you need first to be able to enter complex numbers 
into Mathcad. (This requires a little care, since Mathcad needs to know 
when the symbol i is being used to mean \/—1, rather than in some other 
way, such as, say, part of a variable name.) For example, to enter 3 + 42, 
you can click on the buttons on the ‘Calculator’ toolbar, including the 2 
button, or just type 3+4i. (You do not enter a multiplication between the 
4 and the i.) However, when 7 is entered alone in Mathcad, it must be 
entered as li. This extra ‘1’ in front of the 2 is entered for you 
automatically if you click on the i button, but when using the keyboard, 
you must type 1i. For example, to enter 2 — 7 via the keyboard, you type 
2-11. (The ‘1’ is visible only when entering or editing the complex 
number. It disappears once the number is complete.) 


Activity 6.1 Entering complex numbers in Mathcad 


(a) By hand, simplify each of the following complex numbers. 
(i) (8 + 47)(3 — 42) (ii) 2(3 — 42) 

(b) Create a new (Normal) worksheet. Then enter each of the complex 
products in part (a), and enter = to evaluate them. 


Comment 


(a) These products are 
(i) 2a: (ii) 4+ 32. 

(b) Mathcad gives the above values. If you use keyboard entry, and type i 
or 1*i to enter 7 in (ii), then Mathcad gives an error, treating 7 as an 
undefined variable, not as \/—1. The key sequence required here is 
1i* (3-41) . (Whatever entry method you use for (ii), you must enter 
the multiplication between the 2 and the left bracket — Mathcad will 
not insert it if you omit to do so.) 


The next activity provides further practice in entering complex numbers. 


This symbolic keyword was 
introduced in MST121 
Chapter A3; see also A Guide 
to Mathcad. 


Remember, to enter a + bi, 
use the buttons on the 
‘Calculator’ toolbar or type 
at+b*1i. (Note that, when 
variables are used to 
construct a complex number, 
a multiplication * is required 
before the 71, which must be 
typed as 1i.) 


Remember that names are 
case sensitive in Mathcad. 
You must use capital R and I 
here. 


CHAPTER D/ 


Activity 6.2 Multiplying complex numbers 


(a) By hand, evaluate (1 — 27)(3 + 47). 


(b) On a new (Normal) worksheet, enter the expressions shown below. For 
each of zw and (ac — bd) + (ad + bc)i, click on the expression and then 
enter = to evaluate it. Experiment with other values of a, b, c and d, if 


you wish. 
a= 1 bez c:=3 d:=4 
Zi=at bi weet di 
Zu 
(ac-bd)+f(ad+ boy 
Comment 


(a) This product is 11 — 27, as Mathcad confirms. 


(b) The two expressions give the same value whatever values are specified 
for a, 6, c and d, as one would expect! 


Leave this worksheet open; it will be used again for the next activity. 


Mathcad can evaluate the real and imaginary parts of a complex number, 
complex conjugates and the modulus of a complex number. It uses the 
usual notations: Re(z) for the real part of z, Im(z) for the imaginary part 
of z, Z for the complex conjugate of z, and |z| for the modulus of z. To 
enter the first two, one just types the expressions Re(z) or Im(z). To 
enter Z, type z[Shift]2. To enter |z|, type z, then click on the |z| button 
on the ‘Calculator’ toolbar or use the keyboard alternative [Shift]\. 


Activity 6.3 Complex arithmetic 


(a) With the complex numbers z = 1 — 27 and w = 3+ 4i entered on a 
worksheet, as you did for the previous activity, use Mathcad to 
calculate each of (i)—(ix) below. 


(i) z+w (ii) w+z (iii) Re(z) (iv) Im(w) (vi 2 
(vi) |w| (vii) 1/w (viii) w/z (ix) wz /iz/* 
(b) Try varying a, b, c and d, and check that (viii) and (ix) in part (a) give 
the same answers each time. 
Comment 
(a) Mathcad gives the following answers. 
(i) 442i (i) 442i (ii) 1) (iv) 4 (W) 142) (i) 5 
(vii) 0.12 — 0.162 (viii) —1+ 2% (ix) —1+ 2% 
(b) We have z x Z = |z|? (see Exercise 3.2(b)), and so 
wiz = we /2e = we] {|z). 


Mathcad should therefore give the same values for w/z and for 
wz/|z|*, for any complex numbers w and z (where z # 0). 
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Mathcad notes 


© Re(z) and Im(z) are built-in Mathcad functions. As well as typing 
them directly, you can also enter them by selecting Function... from 
the Insert menu. Choose ‘Complex Numbers’ from the ‘Function 
Category’ box, then either ‘Re’ or ‘Im’ from the ‘Function Name’ box, 
before clicking on ‘Insert’. 


© The double-quote " (obtained by typing [Shift]2) performs several 
roles in Mathcad. When entered after an expression (selected within 
the blue editing lines) it applies the complex conjugate operator to 
that expression. However, when " is entered in an empty space in the 
worksheet (at the red cross cursor) it creates a text region, and when 
entered in an empty placeholder in an expression it creates a text 
string variable. 


We look next at Argand diagrams, and at translating between the polar 
and Cartesian forms of a complex number. Mathcad file 221D1-01 provides 
templates for translating between the two forms. It also shows Mathcad 
plots illustrating each of these forms. 


Activity 6.4 Polar and Cartesian forms 


Open Mathcad file 221D1-01 Cartesian and polar forms, and turn to 
page 2 of the worksheet. Here a and 0 are the real and imaginary parts of 
the complex number z, whose Cartesian form is z =a+ bi. The 
corresponding polar form is (r,@), where r = |z| and 6 is defined by the 
Mathcad expression arg(z), of which more below. In the polar form, @ is 
by default in radians, but @ is also shown in the worksheet in degrees 
(which is usually easier to visualise). The page also shows Argand diagram 
plots illustrating each of the Cartesian and polar forms. 


(a) Working by hand, express each of the complex numbers (i)—(iv) below 
in polar form. 


(i) 2-27 (ii) —5 (iii) 37 (iv) —2 — 22 
(b) Check that Mathcad gives in each case the result for the polar form 


that you calculated in part (a). In what range do the values that 
Mathcad gives for arg(z) lie? 


(c) Use the worksheet to obtain the polar form of 3 — 82. 


(d) Given a complex number (r,@) in polar form, whose Cartesian form is 
a+ bi, what are its real part a and imaginary part 6b? 


(e) Turn to page 3 of the worksheet, which gives a template for converting 
from polar to Cartesian form. Check that the equations used here are 
as you would expect. Use this template to express each of the 
following in Cartesian form. 


(t} 43,9) (i) 12,280 
Comment 


(a) We obtain the following polar forms. 


(i) (2V2,-4m) (ii) (5,0) (ii) (3,57) (iv) (22, — 3) 


Remember that Mathcad 
notes are optional. 


You will not be asked to 
create such Mathcad plots. 
(For your interest only, details 
of how to do so are provided 
on page 4 of the worksheet.) 


Mathcad numbers the angular 
grid lines in a polar plot from 
0° to 360°. 


Solve blocks were used in 
Mathcad file 221B1-03 for 
Chapter B1. The 
Newton—Raphson method was 
applied in Mathcad file 
221C1-03 for Chapter Cl. 


CHAPTER Di 


\ 


(b) Above the Argand diagram plots, Mathcad gives the same values that 
were found in part (a), but expressed as decimals (to 3 decimal 
places), as follows. 


(i)< {258285 D785) 3)-45/8:dd2) 991A) 
(iv) (2.828, —2.356) 


Mathcad also gives exact answers, obtained by using symbolic 
evaluation (—>), further down page 2 of the worksheet. Here J/2 
appears as 21/2, 


The values that Mathcad gives for arg(z) lie in the range (—7, 7]. 
Thus the Mathcad function arg gives the principal value of the 
argument, as defined in Subsection 3.2. (However, when expressed in 
degrees, these values differ by 360° from those shown on the lower half 
of the Argand diagram polar plot.) 

(c) You should obtain the polar form (8.544, —1.212) (to 3 decimal places). 

(d) We have a=rcos@ and b=rsin8. 

(ce) We obtain the following. 
(i) —3 (as we would expect) 


(ii) —1.010 + 1.7262 (to 3 decimal places) 


Now close file 221D1-01. 


6.2 Finding roots with Mathcad 


You have met various ways of using Mathcad to solve equations. Solve 
blocks and the Newton—Raphson method each employ a numerical 
algorithm, and each will find only one solution at a time. To find all the 
roots of a polynomial, both real and complex, a different approach is 
preferable. 


Mathcad will give the two roots of a quadratic polynomial, using the 
symbolic keyword ‘solve’. This works whether the roots are real or 
complex. The roots are obtained using the formula for the solutions of a 
quadratic equation, to give exact rather than approximate solutions. If you 
enter the quadratic with its coefficients given as integers or rationals, then 
Mathcad gives the roots of the quadratic exactly, using square roots if 
necessary. If you enter the coefficients as decimals, then Mathcad returns 
the roots as decimals. Figure 6.1 shows an example. Mathcad initially 
gives these decimal solutions to 20 places, but if you click on the answer 
and enter =, it displays the solutions to the number of decimal places 
specified in ‘Number of decimal places’, on the ‘Number Format’ tab under 
Result... on the Format menu. 


: oe 7 
2x +—x+— solve,x 4 
Z 4 


2 


— A2500000000000000000 -— 6959705453537 5274026 x 1 —0,.625 - 0.6961 
2x + 2.5x+ 1.75 solve,x 4 = 


— 62500000000000000000 + 6959705453537 5274026 x 1 —0.625 + 0.6961 


Figure 6.1 Roots obtained using the symbolic keyword ‘solve’ 
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This distinction, between exact and decimal expressions, may not seem 
very important, but it becomes more significant for a cubic polynomial. 
There is a formula giving the roots of a general cubic (we touched on this 
in Section 1), but it is much more complicated than the formula for the 
roots of a quadratic. Mathcad uses that formula if you use the symbolic 
keyword ‘solve’ to find the roots of a cubic. 


Activity 6.5 Solving a quadratic and a cubic 


Create a new (Normal) worksheet for this activity. 


(a) Use the symbolic keyword ‘solve’ to find the roots of the quadratic 
ae : 
(i) entering the coefficients as rationals; 


(ii) entering the coefficients as decimals. 


With the vertical blue editing line at the right-hand end of the 

expression, click on the ‘solve’ button on the ‘Symbolic’ toolbar, then 
type x into its placeholder and either click elsewhere on the worksheet 
or press [Enter]. Alternatively, just type [Ctrl] [Shift] .solve,x. 


(b) Use the symbolic keyword ‘solve’ to find the roots of the cubic 
=. = 
(i) entering the coefficients as rationals; 


(ii) entering one coefficient as a decimal (for example, enter 2 as 2.0). 


Comment 
(a) You should obtain the results shown in Figure 6.1. 


(b) In (1), you will obtain a large expression that spills off the right-hand 
side of the page and is not easy to read. In (ii), after clicking on the 
answer and entering =, you obtain the roots as —1.521 and 
0.761 + 0.8587 (to 3 decimal places). 


If you click on the large exact expression given in (i) and enter =, 
Mathcad will give the value of that exact solution expressed as a 
decimal to 3 places, and so the same value as in (ii). You may need to 
scroll to the right to see this, though. 


The situation for a quartic polynomial is similar to that for a cubic. There 
is a formula for the roots, and Mathcad employs this when the symbolic 
keyword ‘solve’ is used. The formula again produces very cumbersome 
answers in exact mode, but is usually an effective way of obtaining all four 
roots if the coefficients are entered as decimals. 


For polynomials of order five or higher there is no general formula for the 
roots. In this case, Mathcad may respond in a variety of ways. It is quite 
good at recognising special situations where it is possible to give all the 
roots exactly, but in general this is not possible. Sometimes it gives 
answers with 20 decimal places, even though none of the coefficients is in 
decimal form. 


Mathcad also has an alternative facility, called ‘polyroots’, which 
calculates all the roots of a polynomial, using an iterative procedure. 
Unlike the symbolic keyword ‘solve’, polyroots will not give exact answers, 
but then the exact expressions are often so large as to be unmanageable. 


One advantage of polyroots is 
that it avoids the need to 
enter all the powers of x 
involved into the worksheet. 


To define a matrix or 
subscripted variable, you can 
click on the appropriate 
button on the ‘Matrix’ 
toolbar, or type [Ctrl]m or [ 
(left square bracket ) 
respectively. 
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In most cases, polyroots is the most efficient way of finding all the roots of 
a polynomial. To use polyroots, we need first to enter the coefficients of 
the polynomial into Mathcad, as a vector. To do this, the coefficients can 
be entered either as a matrix with one column, or as subscripted variables. 
For example, to find all the roots of 


a” + 9e* = Fe" a + Se + 1, 
we could use either of the approaches illustrated in Figure 6.2. 


Notice that ag (at the top of the matrix) or bo gives the constant 
coefficient, which is 1 in this example; a, (second down in the matrix) or b; 
gives the coefficient of x (here 2), and so on. Entering =, after clicking on 
either of the expressions polyroots(a) or polyroots(b), gives the roots. 


2 
an by = 1 b, =2 b, = -1 
a= 
—2 
, b,=-2 by=3 bol 
1 
polyrootst a) polyroots(b) 


Figure 6.2 Mathcad screens to find all the roots of x° + 32+ — 22° — a? +2r+1 


Activity 6.6 Finding the roots of a quintic 


Create a new (Normal) worksheet for this activity (or continue with your 
worksheet for Activity 6.5). 


(a) Use polyroots to find all the roots of the quintic polynomial 
g° + 3a* — 20° — o* + Qe + A 
(b) What is the response from Mathcad if you enter the polynomial in 
part (a) and then apply the symbolic keyword ‘solve’, 
(i) seeking exact solutions; 
(ii) seeking approximate solutions, by putting a decimal point into one 
of the coefficients? 
Comment 


(a) Figure 6.2 shows the two approaches that can be used here. We obtain 
the following five roots: 


—3.454; —0.541+0.1617; 0.768 + 0.5662. 


(b) (i) If all the coefficients are entered as integers, then Mathcad’s 
response is to give the roots correct to 20 decimal places. 


(ii) If one or more coefficients is entered as a decimal, then Mathcad 
gives the identical output to that described in (i). 


SECTION 6 Complex numbers and Mathcad 


Activity 6.7 Finding roots 


In parts (a) and (b) below we give two results obtained by hand in 
Section 4 of the main text. Use Mathcad to confirm each of these results. 


(a) In Activity 4.4, you found that 8 has the three cube roots 
2. -1+iv3. 


(b) In Activity 4.7, you found that 42 has the four fourth roots 
1.307 + 0.5417; —0.541+4 1.3072; -—1.307—0.5411; 0.541 — 1.3072. 


Comment 


(a) The cube roots of 8 are the roots of the polynomial z’ — 8. The roots 
have been given exactly, and to obtain exact roots from Mathcad, we 
should use the symbolic keyword ‘solve’. Since we have a cubic 
polynomial here, this should give all the roots, and it does. 
(Remember to enter the coefficients as integers.) 


o 


The fourth roots of 47 are the roots of x* — 47. Use of polyroots gives 

the roots in the form above. (Enter the coefficients aj = —41, If you use the same worksheet 
a, = a, = a3 = O and a, = 1.) The symbolic keyword ‘solve’ displays here as for Activity 6.6, and 
the solutions in an alternative format, which gives the same values to enter the coefficients using 

3 decimal places after entering =. If either coefficient is entered as a separate subscripted 


decimal, then ‘solve’ finds only one of the four solutions. variables, then pole Seay 
inherit the coefficient as = 1 


Mathcad notes from that work, and 
polyroots(a) may not give the 


Care is needed in Mathcad when redefining subscripted variables. Suppose 
answer you require! To avoid 


that you define six variables, ao, a,,...,@5, but then later in the same this problem, use a different 
worksheet, you redefine only five of them, ao,@1,...,a4. The sixth value as name for the subscripted 

is still defined, as it is ‘inherited’ from the first definition. (Note that while variables in this activity — see 
this problem may occur when defining separate subscripted variables, it the Mathcad notes for an 


does not occur if the definitions are made by defining a matrix a with one __ explanation. 
column. ) 


You saw how to calculate the nth roots of unity in Section 4. There, we 
looked at some particular values of n. We can use the same approach to 
find a general formula for roots of unity. The nth roots of unity satisfy 
z” = 1. If we express z in polar form, as z = (r,@) and 1 in polar form, as 
(1,0), this equation gives 


A S60 = ir” , ne). 


For this to hold, we need r = 1, and né to be a multiple of 27. Thus, in 
polar form, the nth roots of unity are (1,2ka/n), where k takes the 
values 0,1,...,2—1. In Cartesian form, they are 


2k 2k 
cos (=) +27sin (=), or t= 0, 1,...., 7 — 1. 
n n 


The next activity features a Mathcad file to find nth roots of unity. 


if 


This subsection will not be 
assessed. 


iZ 


CHAPTER Di 


Activity 6.8 Roots of unity 


Open Mathcad file 221D1-02 Roots of unity, and turn to page 2 of the 
worksheet. This page is set up to calculate the nth roots of unity (for a 
given value of n), both using the formula above and using polyroots, to 
enable you to compare the results. The coefficients aop,a ,,...,@,_1 are 
defined by means of the Mathcad function if(condition, p, q). 


(a) Use this page to find the ninth roots of unity. 


(b) Examine the Argand diagram plot of the nth roots of unity, on page 3 
of the worksheet, for n equal to each of 5, 10, 25 and 50, and note any 
patterns that you observe. (You will need to set the value of n on 
page 2, in each case.) | 


Comment 


(a) You need to set n equal to 9. The formula and polyroots give the same 
roots (though in different orders). The ninth roots of unity are 


1, 0.766 + 0.6432, 0.174 + 0.9852, —0.5 + 0.8662, —0.940 + 0.3422. 


(b) The roots of unity are always evenly spaced around the unit circle, as 
mentioned in Section 4. For larger values of n, the plot of 
corresponding points looks increasingly like a circle. (You may have 
noted other features. ) 


Mathcad notes 


The tables of roots on page 2 of the worksheet have been formatted to 
display them in the same style. By default, entering ‘polyroots(a) =’ 
displays nine values or less within round brackets (in matrix form). To 
force Mathcad to display a table here, for any number of values, we have 
used the Format menu, Result..., and on the ‘Display Options’ tab, 
changed the ‘Matrix display style’ from ‘Automatic’ to ‘Table’. Further 
changes have been made by clicking on the ‘polyroots’ table with the right 
mouse button, and choosing first ‘Properties...’, then ‘Alignment’ from the 
resulting mini-menu. In the table properties, column/row labels are not 
shown, and the alignment (of the expression ‘polyroots(a) =’ to the table 
values) is set to ‘Top’. This alignment has also been set for the other two 
tables on the page, ‘k =’ and ‘z, =’. 


Now close file 221D1-02. 


6.3 Complex exponentials and geometry 


Complex exponentials are entered into Mathcad in just the same way as 
real exponentials, for example, by typing e/z (or exp(z)). 


Activity 6.9 Complex exponentials (Optional) 


(a) In a new worksheet, set z = 2 + 32, then evaluate e’. 


(b) With z = 2+ 32, what are |e*| and arg(e*)? Use Mathcad to verify 
your answers. 
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Comment 


(a) We obtain e?** = —7.315 + 1.0437. (You found this result by hand in 
Activity 5.2.) 


(b) We have e* = e?7* = e?(cos(3) + isin(3)). So 
tae are’ 7=9. 
Mathcad confirms these results. We obtain |e*| = 7.389, and this 


equals e? to 3 decimal places. Symbolic evaluation (—>) of |e*| gives 
exp(2) directly. 


We now use Mathcad to investigate sequences generated by the recurrence 
system 


Co=1, Cri =ke, (n=0,1,2,...), (6.1) 


where k is a fixed complex number. (We considered such sequences in 
Subsection 5.2 of the main text.) 


Activity 6.10 Iterations with |k|=1 (Optional) 


Open Mathcad file 221D1-03 Complex exponentials and turn to 
page 2 of the worksheet. Here we consider the recurrence system (6.1) with 
k of the form e’’, so that |k| = 1. 


(a) Consider first the sequence generated by the recurrence system (6.1) 
with k= e'/*, 
(i) This sequence will eventually repeat itself. How do we know this? 


(ii) How many iterations are needed before the sequence starts to 
repeat? Set N on page 2 of the worksheet to the smallest value that 
will plot all points of the sequence. Check that larger values of N plot 
nothing new. 


(b) For what values of 6 in k = e’® does the sequence generated by the 
recurrence system (6.1) eventually repeat? 


(c) Decide on a value of 6 for which the plot should look close to a circle. 
Choose a value of 6 for which the sequence eventually repeats, and 
choose N large enough to show all the points of the sequence. Enter 
these values on page 2 of the worksheet, and check that they have the 
effect that you expected. 


(d) Decide on a value of 6 for which the sequence does not eventually 
repeat. Examine the plot on page 2 for your chosen value of 6. 


Comment 


(a) (i) You saw in Activity 5.5 of the main text that the sequence 
generated by the recurrence system (6.1) eventually repeats if and only 
if k is a root of unity. Now 


(ein /S)? = eliniei =e" = cos 2n) + tsin{ 2a) = 1, 


SO eit /6 


(ii) We have cy. = k'? = 1 = cp, and after this the terms of the 
sequence repeat. So with N = 11 on page 2 of the worksheet, we see 
all points of the sequence plotted. With N = 12, the line joining c,; 
and cij2 completes a regular 12-sided polygon, and the Mathcad plot 
looks the same for all higher values of NV. 


is a twelfth root of unity. 


a3 


You should still be working 
with Mathcad file 221D1-03, 
on page 2 of the worksheet. 
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(b) 


The sequence eventually repeats if k is a root of unity. Now k = e” is 
a root of unity if k” = 1 for some integer m; that is, if we can find an 


integer m such that 
Cad wee ee ee 


This is the case if m@ is equal to 27, or to an integer multiple of 27. 
Thus the sequence generated by (6.1) with k = e’” eventually repeats if 


6 = 2pr/m, 
where p and m ( 0) are integers. 


The plot will eventually repeat if we choose k = e”” with 0 = 2pm/m. 
It will look close to a circle if we choose m reasonably large, say 

m = 100, and for simplicity we may as well take p = 1. So take k equal 
to, say, e’"/°° and N = 99. (There are many other possible choices of k 
and JN.) 


Any value of k that is not a root of unity will produce a sequence that 
does not repeat. For example, for k = e’ with 6= 1, and N =7, we 
obtain the plot shown in Figure 6.3(a). You can see that k’ is not 
equal to 1, and the plot does not return to its starting point. Nor does 
it return to 1 later. For example, the plot with N = 13 is shown in 
Figure 6.3(b). With larger values of N, as illustrated in Figure 6.3(c) 
with N = 100, it is much harder to make out what is happening, but 
the iteration never returns exactly to its starting point of cp = 1. 


Figure 6.3 Plots of the sequence generated by recurrence system (6.1) with 
k = e’, showing the sequence up to k% with N equal to: (a) 7 (b) 13 (c) 100 


Mathcad notes 


To enter the complex exponential re’? in Mathcad, you can use the buttons 
on the ‘Calculator’ and ‘Greek’ toolbars, or type r*eA1i*q[Ctrll]lg. (You 
must type 1i and include the multiplication * between the 7 and the 0.) 


Activity 6.11 The twentieth roots of unity (Optional) 


(a) 
b) 


er, 


For what values of @ is e*® a twentieth root of unity? 


Examine the sequences generated by the recurrence system (6.1) with 
k = e'P™/10 where p is equal to each of 


Linky 8825) 6% TA 


Note any points that you observe about the plots you obtain. 


SECTION 6 Complex numbers and Mathcad 


Comment 


(a) We have (e’%)?° = 1 if e°° = 1. This is the case if 200 is an integer 
multiple of 27; that is, if 9 = px/10, where p is an integer. (We obtain 
the complete set of twentieth roots of unity if p takes the values 
0,1,2,...,19; other values of p just repeat one of these.) 


(b) Plots for p = 1, 2 and 3 are shown in Figure 6.4 (for N = 20 
iterations). For p = 1 we see the twentieth roots of unity uniformly 
spaced around a unit circle, as we would expect. For p = 2, we obtain 
fewer points, 10 rather than 20. This is because 2(7/10) = 7/5, and 
(e7/5)!0 — e?™ — 1. Hence e?*/1° is a tenth root of unity (as well as 
being a twentieth root). So the plot in this case repeats itself after 
10 iterations, and not all the twentieth roots of unity are visited. With 
p = 3, all the twentieth roots of unity are again visited by the plot. 
The order in which they are visited is different, though. With p = 1, 
the roots are visited in order as we go round the circle, that is, with m 
equal to 0,1,2,...,19 in e*”"/7°, With p = 3, the roots are visited in 
the order 


Oe a 0 9 119, 2, GS 8 11, 14, 17. 


Figure 6.4 Plots of sequences generated by recurrence system (6.1) with 
k = e?7/10 for p equal to: (a) 1 (b)2 (c)3 


With p = 4, we see only five roots. Here 4(7/10) = 27/5, and e”'”’° is 
a fifth root of unity. 


With p =5, we see only four roots, since 5(7/10) = 7/2 and e’””? is a 
fourth root of unity. 


With p = 6, we see 10 roots, since 6(7/10) = 37/5, and e**”/° is a 
tenth root of unity. Here the roots are visited in the order (labelling 
round the circle): 0, 3, 6, 9, 2, 5, 8, 1, 4, 7. 


With p = 7, all 20 roots are visited. Notice that 77/10 does not 
simplify at all, since 7 and 10 have no common factors. So here e 
is a twentieth root of unity, but not an nth root for any smaller value 
of n. 


Tin /10 


With p = 17, we see the same plot as for p = 3. Although you cannot 
see it from the plot, in this case the roots are actually visited in the 
reverse of the order for p = 3, that is, 0,17,14,11,...,3. 


With p = 18, we see the same plot as for p = 2; and for p= 19, we 
obtain the same plot as for p = 1. 


If |k| is not equal to 1, then the recurrence system (6.1) produces a spiral, 
as you saw in Subsection 5.2. 


15 


You should still be working 
with Mathcad file 221D1-03, 
on page 2 of the worksheet. 


You should still be working 
with Mathcad file 221D1-03. 
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Activity 6.12 Spirals (Optional) 


Consider now sequences generated by the recurrence system (6.1) for a 
general complex number k, with polar form (r,@) where r F 1. 


(a) Examine the plot of the sequence with the following values of 


k = tee). 
(i) (1.1,7/6) (ii) (V2,/2) (ii) (1/2, 1/2) 
(iv) (V2, —m/2) 


In each case try two or three values for N. You will need to adjust the 
value of the graph scale variable s appropriately for each value of k. 


(b) Choose your own values of r and @ in k = (r,@). How does the spiral 
vary? 
Comment 


(a) You saw plots of the spirals in these cases in Figures 5.2 and 5.3 of the 
main text. 


(b) With r > 1, the plot spirals outwards, and the larger r is, the more 
rapidly the spiral moves out. If r < 1, the spiral moves inwards. If 
0 <6@< 7 the spiral goes round anticlockwise, while if —7 < 9 < 0 it 
goes round clockwise. 


Finally, we look at plots of the complex-valued function f(t) = r(t)e”. 
With r(t) =a‘, this plot-gives a continuous spiral. 


Activity 6.13 Continuous spirals (Optional) 


Turn to page 3 of the worksheet, which shows a plot of the complex-valued 

function f(t) = r(t)e” (t > 0). 

(a) Examine the plot with r(t) = a’, where a is equal to each of the 
following. 


(i) 12 jay 12 (iii) 0.9 (iv) 0.8 (v) 1 
How does the plot vary? (You will need to adjust the value of the 
graph scale variable s in order to see the plots for (iii) and (iv) well.) 


(b) Examine the plot with r(t) equal to each of the following. 
(i}-Z (ii) 1+0.2cost (i) Int GS 1) 
How does the plot vary? (Again, adjust s appropriately. ) 


Comment 


(a) With a > 1, the plot is an increasing spiral. The rate of increase is 
faster if a is larger. With a < 1, the plot is a decreasing spiral. The 
rate of decrease is faster if a is smaller. With a = 1, the plot is a circle. 


(i) With r(t) = t, we again see an increasing spiral. 


S 
4 


(ii) In this case the plot is a closed curve. 


(iii) Here you need to change the range for t, setting T1 := 1, because 
the domain given is t > 1. We see an increasing spiral, but the rate of 
increase becomes slower and slower. (This is particularly clear if you 

modity the range for t, to give a larger upper limit, such as T2 := 40.) 


Now close file 221D1-03. 


Chapter D2, Section 5 
Number theory and Mathcad 


In this section, you will have an opportunity to develop your 
understanding of the concepts in the main text, to check some of your 
earlier calculations, and also to try out the methods with much larger 
numbers than your calculator can handle. 


Mathcad can perform some arithmetic operations with very large integers 
by using symbolic evaluation, but we do not describe these capabilities 
until the end of the section. 


5.1 The Division Algorithm 


The Division Algorithm asserts that if a and n are integers, with n See Theorem 1.1 of the main 
positive, then there are unique integers g and r such that text. 
a=qn+r, withO<r<n. (5.1) 


The number gq is called the quotient and r is called the remainder of a, on 
division by n. 


Equation (5.1) leads to formulas for q and r: 


q = floor (=) and r=a-—qn. As you saw in Subsection 1.1 
‘i of the main text, floor(z), 
The first activity in this section asks you to try out these two formulas in also denoted by [z], is the 
Mathcad. greatest integer less than or 
equal to x. 


Activity 5.1 Quotients and remainders 


(a) By hand, find the quotient and remainder of 
(i) 32 on division by 6; (ii) —32 on division by 6. 


(b) Create a new (Normal) worksheet, and enter the following expressions. 


aie 32 i= 6 


a 
= floor] — r=a-qn 
é (=) ‘ Make sure that you place the 
formula for r after that for q 


Then enter q= and r= to evaluate these quantities. Change the value in the worksheet. 


of a to —32, and check that Mathcad gives the same answers that you 
obtained in part (a). 

Comment 

(a) ti) Ga, f= 2 (ii) q=-6,r=4 

(b) Mathcad gives the same values. 


It is a good idea to enter some text to structure the worksheet you are 
creating in this and subsequent activities. For example, a title and some 
headings and comments might be useful. 


ly 
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Activity 5.1 shows that Mathcad can easily be used to find quotients and 
remainders. To investigate the behaviour of these quotients and 
remainders, it is convenient to introduce a pair of functions defined as 
follows. 


quotia,n) = foo * rem an) = a- nquotta,n) 
fi 


Activity 5.2 Varying a and keeping n fixed 


(a) Enter the functions above in your Mathcad worksheet, and check that 
quot(—32,6) = —6 and rem(—32, 6) = 4. 


(b) Describe how quot(a,n) and rem(a,n) behave when n is kept fixed 
and a is allowed to increase. It may help to consider an example, such 
as n:= 3 and a := —10,—9..10, and to create tables or graphs for 
quot(a,n) and rem(a,n). (You may need to format graphs 
appropriately to show the behaviour clearly.) 

Comment 

(b) As a increases 


© quot(a,n) increases by 1 whenever a reaches a multiple of n, but is 
otherwise equal to its previous value; 


© rem(a,n) cycles repeatedly through the numbers 0,1,2,...,n—1. 
These statements are illustrated by the graphs in Figure 5.1. 


5 


quotta,n) oO 
ae 4 


Figure 5.1 Graphs of the functions (a) quot, (b) rem 


You may have wondered whether Mathcad has its own built-in functions 
for calculating quotients and remainders. In fact Mathcad has a function 
called mod, which can be used for finding remainders. This function has 
two arguments: it can be obtained by typing mod(a,n). 


Activity 5.3 The function mod 


(a) Enter the following expressions in your worksheet, and then evaluate 
them by entering =. 


(i) mod(32, 6) (ii) mod(—32, 6) 


(b) Vary the numbers used as arguments in part (a). On the basis of the 
observed values, attempt to define the function mod. 


(c) What happens when you try to evaluate mod(a*,n), where a = 3, 
k = YOOO and # = 137 
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Comment 
(a) (i) mod(32,6) = 2 (ii) mod(—32,6) = —2 
(b) Several illuminating examples are: 
mod(0,6) =0, mod(32,—6)=2, mod(—32, —6) = —2. 


It seems from these examples that Mathcad evaluates mod(a, 7), 


where a and n are integers with n # 0, by first finding the remainder Mathcad will also return a 
when |a| is divided by |n|, and then giving this remainder the same value for mod(a,n) when a 
sign as a. and n are not integers, but 


= ie ; that d t here. 
(c) An error message is displayed: ‘Found a number with a magnitude ce ne er ee 


greater than 10307 while trying to evaluate this expression.’. Finding 
such remainders, when a* is large, is addressed in Activities 5.4 
and 5.6. 


Save the worksheet that you have created, if you wish. Then close the file. 
An alternative definition of the Mathcad function mod is as follows: 
mod(a,n) is the unique integer r satisfying 
o |r| < (nj; 
© az=qn+r, for some integer q; 


© r has the same sign as a. 


This makes it clear that mod(a,n) is a remainder of a on division by n, 
but for a < 0 it is not the remainder which is useful in number theory. 
However, we can use mod(a,n) to find such remainders when a and n are 
both positive. 


Warning: As you will see in Subsection 5.4, when evaluated symbolically 
(with —) the Mathcad function mod does behave in the way required for 
number theory, always giving non-negative remainders! 


5.2 Remainders of powers 


In Section 1 of the main text two methods, here called algorithms, were 
described for finding the remainder of a* on division by n, where a, k and 
n are positive integers. Both these algorithms — repeated multiplication 
and repeated squaring — are implemented in Mathcad file 221D2-01. 


Repeated multiplication 


To find the remainder of a* on division by n, we calculate successively the 


remainders of a°,a', a”, a®,... on division by n, using a recurrence relation The sequence of remainders 
which expresses each remainder in terms of the previous one. always shows a repeating 
pattern. 


If we denote by r; the remainder of a’ on division by n, then 
r, =a‘ (mod n), 

SO 
r, =ax.a'.' =ar;_, (mod n). 


Also, ro = a° = 1 (mod n). Thus, there is a simple recurrence system for 
calculating the remainders 79,171,1T2,.--,Tk- 


19 


You should still be working 
with Mathcad file 221D2-01, 
on page 2 of the worksheet. 


27” = 1024 = 3 x 341+1 


20 


CHAPTER D2 


Activity 5.4 Repeated Multiplication Algorithm 


Open Mathcad file 221D2-01 Remainders of powers, and turn to 

page 2 of the worksheet. This implements the Repeated Multiplication 
Algorithm by means of a Mathcad program. The page is set up with 

a = 3, k= 1000 and w=13, and shows.that the remainder of 2° = 3!" on 
division by 13 is r, = 3. The table of remainders illustrates the repeating 
pattern of values within the sequence 79,11,12,.--5Tr- 


(a) Use the worksheet to find the remainder on division of 5?°°° by 21. 


(b) Experiment with various choices of a and n (keeping k fixed at 2000), 
to find repeating sequences of remainders with the following forms. 


iy 1,6 £97649 By x. (ae 34,1, 7,8, 6, x. 


Comment 


(a) With a = 5, k = 2000 and n = 21, the required remainder is rgo99 = 4. 
Note that the table shows the sequence 


1, 5, 4, 20, 46, 17.1, 5,4. eee. 


with a pattern which repeats after 6 terms. If you scroll across the 
table to 2 = 2000, you will see that the remainder rz999 = 4 appears 
there as well. 


A ‘by hand’ check is as follows. Since 5° = 1 (mod 21), we obtain 
5S ee = (5) x 5° SB = 25 = 4 (mod 21). 
(b) (i) a=4, n = 12 gives remainders 1,4,4,... 
(ii) a = 3, n = 12 gives remainders 1,3,9,3,9,... 
(iii) a = 2, n = 14 gives remainders 1, 2,4, 8, 2,4,8,... 


These answers are not unique. 


The next activity revisits a result from Section 3 of the main text. 


Activity 5.5 Experiment with primes 


(a) Use the worksheet again to experiment with various prime numbers 
for n, and positive integers a. (Keep k = 1000.) Notice from the table 
that 1 usually appears in the sequence of remainders, and try to 
explain this. 


(b) Try n = 341 (which is not prime) and a = 2, and check that in this 
case also, 1 appears in the sequence of remainders. 


Comment 


(a) If n is a prime number p and a is not a multiple of p, then we expect 1 
to appear in the sequence of remainders, by Fermat’s Little Theorem 
(Theorem 3.2 in the main text): 


Let p be a prime number, and let a be a positive integer which is 
not a multiple of p. Then a?~! = 1 (mod p). 


(b) You should have found from the table that 2'° = 1 (mod 341), which 
can also be checked by hand. ‘This shows that 1 may appear in the 
sequence of remainders, even when n is not a prime number. 
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Repeated squaring 


A more efficient algorithm for finding the remainder of a* on division by n 
involves the steps set out below. (Some of the details are illustrated in the 
margin for the case a* = 147", n = 55, discussed in Subsection 1.3 of the 
main text.) 


© Represent k as a sum of powers of 2: 


k= +a xt 2 te XS se xo”, (5.2) 


where €);¢),...,¢, are in. Z, and c,, = 1. 


© 


Find the remainders S89, 81, 82,..-,5m On division by n of 
a',a?,a*,...,a” , using repeated squaring: 


8; = s7_, (mod n). 


© 


Find the remainder on division of a* by n: 
x ta ¥e 


52 ae ee ese (aed) 


ak Se yg? cade (a?) x (a*)° ene 


co 


= C1 
= Sy 


xX sf 
using repeated multiplication modulo n. 


This algorithm is implemented on the next page of the Mathcad 
worksheet. It is a little complicated, particularly the first step which finds 
the representation of k as a sum of powers of 2, the so-called binary 
representation of k. We shall explain how the implementation of the 
algorithm works after you have used it in the next activity. 


Activity 5.6 Repeated Squaring Algorithm 


Turn to page 3 of the worksheet. This implements the Repeated Squaring 
Algorithm by means of a Mathcad program. The page is set up with 

a = 14, k = 27 and n = 55, and shows that the remainder of a* = 14?" on 
division by 55 is p,, = 9. The table displays intermediate values generated 
by the algorithm. 


Use the worksheet to find each of the following: 
(a) the remainder of 2'°°°° on division by 10001; 


(b) the remainder of 2°°' and of 3°°* on division by 561. 


Comment 
(a) Here you need to enter a = 2, k = 10000 and n = 10001. The 


algorithm gives the remainder as p,, = 4674, so 
210000 = 4674 (mod 10001). 


(This result shows that 10001 is not a prime number, because if it 
were then the remainder would have been 1, by Fermat’s Little 
Theorem. As you saw in Subsection 3.3 of the main text, we have 
10601 = 7a x a) 


Here the remainders p,,, are 2 and 3, respectively, so 
2°61 = 2 (mod 561) and 3°°' = 3 (mod 561). 


3 


(This result suggests that a°°* = a (mod 561), for every positive 
integer a, even though 561 is not prime, and this is true. A proof of 
this fact can be based on the congruences 


(mod 3), a’ =1 (mod 11), 
which hold if a is coprime with 561.) 


a = 1 Geod 17), 


= 


27 = 14248416 
= 141x230 x? 
3 ea) eo 


For a = 14 we have, for 
example, s; = 31 and so 


$2 = 31° = 26 (mod 55). 
14! x 142 x 148 x 14/6 


14! x 31! x 16! x 36 
(mod 55) 


1427 


You should still be working 
with Mathcad file 221D2-01. 


The role of the sequence 
0,P1,+++;Pm, 1s described 
after this activity. 


Si=3x<iix« 17 
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The rest of this subsection 
will not be assessed. 


You are not expected to be 
able to construct such an 
algorithm, but if you have 
time then try to see how it 
works. ~ 


ZZ 


CHAPTER D2 


We now give an explanation of the program for the Repeated Squaring 
Algorithm, given on page 3 of the worksheet. 


The input values a, k and n are declared for use in the program. The 
method of calculating the numbers co, c;,..., Cm, in the binary 
representation (5.2) of k, differs from that used in Subsection 1.3 (it avoids 
having to find the largest power of 2 which is less than k). First note that 
qo = k and cy = mod(k, 2). 


The numbers c;,C2,...,Cm are found by calculating the repeated quotients 
on division by 2: 


gq, = floor ($40) =C, +g X 26 4+-+++¢,2™ +", soc, = mod(q,2); 
qo = Roor{ egy) = + -"4-G,2”*, sO Cg = mod(qo, 2); 
Gm, = floor ($¢n=1) = ms MG mvdal( gq... 2). 


The counter 7 keeps track of how many iterations of this type take place. 
The final iteration gives g; = 1, and the variable m is defined as the 
corresponding value of 7, so m denotes the number of iterations that have 
occurred. 


These iterations take place within a ‘while loop’, which is the block of 
program steps with a solid vertical line to the left and headed by the line 
‘while gq; > 1’. Also within this loop, Mathcad calculates by repeated 
squaring the remainders on division by n of a?,a*,...,a” , which are 
denoted respectively by s,,$2,..., 5m. (The remainder of a on division by 
n is given, before the while loop starts, by sy) = mod(a, 7).) 


The values 89, 51,...,5m are used to calculate the sequence of products 
ee se le c eee ¢ cn 
PPE sy PRA Sy a ss Pr eX Se 


modulo n, again by iteration. The final product p,, is the required 
remainder. 


The program also outputs the values of m and ¢;, s;, p; (¢ = 0,1,...,m), 
so details of the calculation can be displayed in a table. 


Now close file 221D2-01. 


5.3 Arithmetic in Z,, 


Mathcad’s mod function makes it easy to obtain addition and 
multiplication tables for Z, = {0,1,2,...,n—1} since, for example, 


a+,b=mod(a+b,n), whena,be€ Z,. 
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Activity 5.7 Addition and multiplication in Z,, 


Open Mathcad file 221D2-02 Modular arithmetic, and turn to page 2 
of the worksheet. 


(a) Remind yourself of the overall pattern in addition tables for Z,,. Try 
changing n to 9, 10 and 12, say. 


(b) Change addition to multiplication in the definition of p,,, and use the The multiplication will be 
table for multiplication in Z,, to find the multiplicative inverse of 5 displayed as a x 6, since ‘x’ is 
ie oa the default (Math menu, 

Options... , ‘Display’) for 

viewing multiplications in this 

worksheet. 


(c) Use the multiplication table for Z¢ to find the multiplicative inverses 
of 7 and 9 in LZ 96. 


(d) Which rows of the multiplication table for Zo, include 1 and therefore 
all of Zog? What do you notice about each of these row numbers and 
the number 26? Relate your observations to Theorem 3.1 of the main 
text. 


Comment 
(a) The addition table for Z,, has the ‘constant diagonal’ pattern. 
(b) The multiplicative inverse of 5 in Z,, is 9. 


(c) The multiplicative inverse of 7 in Z6 is 15, and the multiplicative 
inverse of 9 in Zo¢ Is 3. 


(d) The rows of the multiplication table for Zyg which include all of Zoe 
are: 1, 3, 5, 7, 9, 11, 15, 17, 19, 21, 23, 25. These are the only integers 
in Z», which are coprime with 26. These observations verify 
Theorem 3.1 of the main text for the case n = 26. 


Finding multiplicative inverses in sets Z, is important, for example in 
Section 4 on cryptography. If n is not too large then multiplicative inverses 
can be read off from the appropriate multiplication table, but when n is 
larger a method based on Euclid’s Algorithm is available. 


Suppose that we wish to find the multiplicative inverse of a in Z,,, where a 
and n are coprime. Euclid’s Algorithm involves 


© applying the Division Algorithm first to n and a, and then repeatedly 
to pairs of remainders until the remainder 0 occurs: 


t= ge +7) 
Qa = get1 +72 
Ty = Q3T2 + T3 


(Buck) 
Pm—-2 = ImlTm-1 aie rm 
es ee Qm+ilm 
(since a and n are coprime, we must have r,, = 1); 
© eliminating the remainders r,,7r2,...,%m—1 from all but the last of the 


above equations, to obtain an equation of the form 
ab=kn+1, with bin Z,. 
Then 6 is the required multiplicative inverse of a in Z,,. 


The algorithm for multiplicative inverses is the subject of the next activity. 


Zo 


You should still be working 
with Mathcad file 221D2-02. 


The rest of this subsection 
will not be assessed. 


Again, you are not expected 
to be able to construct such 
algorithms. 
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Activity 5.8 Multiplicative Inverse Algorithm 


Turn to page 3 of the worksheet. This implements the Multiplicative 
Inverse Algorithm by means of a Mathcad program. The page is set up 
with a =7 and n = 26, and shows that the multiplicative inverse of 7 in 
Z55 is b= 15. The tables display intermediate values generated by the 
algorithm. 


Use the worksheet to find the multiplicative inverse of each of the following. 
(a) 8 in L149 (b) 19 in LZ 100 (c) iim) L a9 139 


Comment 


(a) Setting a equal to 8 and n equal to 19 gives the multiplicative inverse 
of 8 in Zy9 to be 12. 


(b) The multiplicative inverse of 19 in Zio9 is 79. 
(c) The multiplicative inverse of 75 in Z49139 is 10 483. 


An alternative method of finding the answer to part (a) is to use the fact 
that 19 is prime, so 


8'* = 1 (mod 19), 


by Fermat’s Little Theorem. Hence the multiplicative inverse of 8 in Zig is 
congruent to 8!” modulo 19, and this can be calculated using, for example, 
the Repeated Squaring Algorithm. A similar approach can be used in 

part (c), because 49 139 is prime. 


We now explain the program for the Multiplicative Inverse Algorithm, 
given on page 3 of the worksheet. The quotients and remainders in Euclid’s 
Algorithm are found by putting ro = a and then using the iteration 


Tr, = mod(n, a), q, = floor (*) . 
a 


tT; = mod (rj_2, Ti-1), qi = floor (=| 


T§=1 


The iteration continues while 7; > 1, after which m is put equal to the final 
value of 2. 


We know that r,, = 1 (provided that a and n are coprime, as assumed). 
The other remainders r; are eliminated from equations (5.3) by 
multiplying by suitable constants c,,C2,...,Cm41, and then adding up all 
m +1 equations. To eliminate the r; (for 7 < m) we require the constants 
to satisfy 

C42 = Qj+1Cj+1 TF Cj: 
Thus we want 

t= —“Gaptps Pie, “tor 7 = Z,... ee I. (SA) 


The remaining equation, of terms obtained by summing m + 1 equations 
as above but not eliminated by equations (5.4), is 


CN + CoG = C1Q14 + (Cm + Cm419m41)lm- 


We know that r,, = 1, so if we choose c,, = 1 and c,,,; = 0, then this 
becomes 


qnt+ca=cjqat+i1;_ that is, 
a(—c1q, + C2) = (—c1)n +1. 
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Defining cp = —qic; + C2 gives acp = 1 (mod n), but it also gives an 
equation that fits the pattern of equations (5.4) for 7 = 0. Thus we have 


the recurrence system 
es oe ee eC 

6 = —-QGastaet Ca, — foe 7 Se 1, = 2g 0. 

This recurrence system is implemented in the Mathcad program. 


Since acy = 1 (mod n), the multiplicative inverse b of a in Z,, is congruent 
to co modulo n; that is, 


b= cy —n floor (©). 
n 


The value of 6 is output by the program. The values of m, r;, q; and c; are 
also output, so the details of the calculation can be displayed in two tables. 


Now close file 221D2-02. 


5.4 Multiple precision arithmetic 


In Section 4 of the main text, you saw how arithmetic with large integers 


is used in cryptography. In particular, you saw the relevance of large prime 


numbers to a method of public key cryptography based on Fermat’s Little 
Theorem. Symbolic evaluation in Mathcad can be used to perform 
arithmetic with large integers, so-called multiple precision arithmetic, and 
also to discover large prime numbers. 


Arithmetic 


For most numerical calculations, Mathcad maintains 15 digits of precision. 
However, for symbolic evaluations, basic arithmetic operations can be 
performed with rational numbers to many hundreds of digits of precision. 


Activity 5.9 Arithmetic and mod with the symbolic processor 


Create a new (Normal) worksheet. Enter each of the following expressions, 
then evaluate it symbolically, either by clicking on the > button on the 
‘Symbolic’ toolbar or by typing [Ctrl]. , the keyboard alternative. 


(a) 111111111111 + 333 333 333 333 
(b) 111111111 111 — 333 333 333 333 


(c) 111111111111 x 333333 333 333 


q) {tutti 114 
333 333 333 333 

fey Fe 

(i) mod (2!0°°, 10001) 


es toa (h) mod(—7, 4) 


Care is needed in Mathcad 
when entering large numbers 
and reading them in answers. 
The digits are displayed in a 
long string — Mathcad does 
not help by grouping them in 
threes, with gaps between, as 
is usual in text. 


To enter the factorial in 
part (f), click on the n! 
button on the ‘Calculator’ 
toolbar, or type ! (given by 
[Shift]1). 


£0 


In this activity, use the 
worksheet that you created 
for Activity 5.9. 


If you are tempted to try 
calculations like these, then 
remember to save your work 
at each stage, because 
Mathcad may be unable to 
cope with huge calculations. 
(Mathcad’s symbolic 
evaluations can handle 
numbers with about 10000 
digits. ) 
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Comment 


The answers provided by Mathcad are as follows. 
a) 444444 444 444 (b) —222 222 222 222 


37 037 037 036 962 962 962 963 (d) + (e) 1125899 906 842 624 


3 


) 
) 15511210 043 330 985 984 000 000 (g) 2 


20 


( 
(c 
(f 
(h) 1 (Check that evaluating mod(—7,4) numerically gives the value —3.) 
(i) 4674 


Factorisation 


A positive integer can be expressed as a product of powers of prime 
numbers, by using the symbolic keyword ‘factor’. 


Activity 5.10 Factorisation 


Enter each of the following integers, and apply the symbolic keyword 
‘factor’ to it. To do this, either click on the ‘factor’ button on the 
‘Symbolic’ toolbar, followed by pressing [Backspace] twice to remove the 
unwanted placeholder, or type [Ctrl] [Shift].factor . 


(a) 10001 (b) 3220422643 = (c) 50! 


Comment 
(a) 10001 = 73x 137 = (b) 3220422643 = 65537 x 49139 


(aera ee 1 1 ee eee 
x 37 x 41 x 43 x 47 


There is much interest in discovering large prime numbers. By applying 
symbolic evaluation in Mathcad to suitable large integers, large primes 
may be found. A useful rule of thumb here is that 


if an integer n has only small prime factors, then some nearby integer 
may well have large ones. 


So it is a good idea to try factorising numbers of the forms 2” + 1, 3" +1, 
etc. For example, applying the symbolic keyword ‘factor’ to 2*° + 1 gives 


2° 4+1= 3° x 11 x 19 x 331 x 18837001; 
hence 18 837 001 is prime. 


The largest known prime number (in 2003) is 2'°*°°°'” — 1, which has more 


than four million digits. It is one of the family of Mersenne prime 
numbers. These are all of the form 2? — 1, where p is a prime number, and 
they are found by applying a special test to numbers of this form. 


In general, it is much harder to factorise large numbers than it is to 
perform basic arithmetic operations with them. Some large numbers, like 
those of the form 10”, factorise very quickly, but most do not. As you saw 
in Section 4 of the main text, this fact makes it possible to base ciphers on 
large numbers which are the product of two large prime numbers. 


Save the worksheet that you have created, if you wish. Then close the file. 


Chapter D3, Section 1 
Symmetry 


1.3 Using symmetries 


In Section 1 of the main text, you saw how the pattern, or structure, of a 
plane set X can be described by finding its set of symmetries S(X). In the 
optional Mathcad files for this chapter you will see how knowing the 
symmetries can be helpful if you wish to plot the plane set using a 
computer package. These files are self-contained and may be worked 
through after reading the following introductory remarks. 


First, in file 221D3-01, we explain how Mathcad can be used to plot a path 
made up of line segments placed end-to-end. Such a path is called a 
piecewise linear path, or PLP for short, and it is defined by prescribing a 
finite sequence of points in the plane, say p(i), 7 = 0,1,...,n and then 
joining p(0) to p(1), p(1) to p(2) and so on, using line segments. 


P(2) p(n) 


Figure 1.20 A piecewise linear path 


In this diagram, the arrows show the direction in which the line segments 
are plotted. Notice that a PLP with n line segments requires n + 1 points, 
or vertices, to determine it. 


One way to prescribe the vertices of a PLP in Mathcad is to take z to be 
the range variable 7 := 0,1..n, and define p(i) to be a function of 7 whose 
values are vectors, that is points in the plane. The corresponding PLP can 
then be plotted using an X-Y Plot, with the components p(i)) and p(i); as 
the arguments on the x-axis and y-axis, respectively, and with the trace 
type set to ‘lines’. 


In file 221D3-02 we use Mathcad to plot a symmetric PLP corresponding 
to a vector-valued function SP(j,7) of two range variables 7 := 0,1..m—1 
and 2 :=0,1..n —1, by going through the range variables in the following 
order. 


(0,0) —— oe ee ee mess ——-+» (0,n-1) 
(1,0) ——-» (1,1) —~» (1,2) —> ——» (l1,n-1) 


This subsection will not be 
assessed. 


af 
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This facility makes it possible to plot plane sets with many symmetries, 
such as the snowflake and rose window shown below. We do this in files 


221D3-02 and 221D3-03. 


(a) 


Figure 1.21 (a) A snowflake (b) A rose window 
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